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Abstract 

Differentiating matrices allow the numerical differentiation of functions 
defined at points of a discrete grid. The present work considers a type of 
differentiating matrix based on local approximation on a sequence of sliding 
subgrids. Previous derivations of this type of matrix have been restricted to 
grids with uniformly spaced points, and the resulting derivative 
approximations have lacked precision, especially at endpoints. The new 
formulation allows grids which have arbitrarily spaced points. It is shown 
that high accuracy can be achieved through use of differentiating matrices on 
non-uniform grids which include "near-boundary 11 points. Use of the 
differentiating matrix as an operator to solve eigenvalue problems involving 
ordinary differential equations is also considered. 

Submitted for publication to the International Journal for Numerical Methods in Engineering. 
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Introduction 


The present work on differentiating matrices has its origins in the study 
of the differential equations associated with rotating beam configurations 
which model the vibrations and aeroelastic stability of rotating structures 
such as helicopter rotor blades and propeller blades. The fourth-order 
boundary value problems associated with rotating beams do not, in general, 
have useful closed-form solutions. Consequently, most theoretical work on 
these problems has been asymptotic or numerical in nature. 

In one approach to numerical solution of these problems, harmonic time 
dependence is assumed to reduce the governing partial differential equation to 
an ordinary differential equation in space which involves an eigenvalue. The 
fundamental derivative representing beam curvature is then taken as a new 
dependent variable, and the eigenvalue problem is rewritten as a second-order 
integro-differential equation (White & Malatino, 1975). This formulation can 
be conveniently expressed using appropriate integral and differential 
operators. In general, boundary operators will also be required. 

The operator equation for the continuous solution may be converted into a 
matrix operator equation for a finite dimensional solution vector by 
evaluating the operator equation at a finite set of discrete grid points which 
span the interval of interest. The key question now becomes the manner in 
which the differential and integral operator matrices are approximated. 

The present work examines one approach which approximates a differential 
operator through use of a differentiating matrix. The specific type of 
differentiating matrix considered uses a sequence of sliding subgrids to 
obtain local approximations at grid points. To maintain simplicity while 
illustrating the basic features of this technique, the discussion here will be 



restricted to considering local approximations obtained through Lagrange 
interpolation. Similarly, the eigenvalue problems considered here are the 
simple historical test problems involving the harmonic equation. Applications 
of the present technique to a variety of rotating beam problems will be given 
in a later paper. 

Differentiating matrices provide a means of numerically differentiating a 
function that is expressed in terms of the values of the function at 
increments of the independent variable, i.e., on a discrete grid with a finite 
number of points. When the values of the function at the grid points are 
arranged as a column vector, premultiplication by a differentiating matrix 
produces a vector containing approximate values of the derivative at the grid 
points. 

An attractive feature of differentiating matrices is that their 
derivation requires only information about the grid points, and no explicit 
information is needed about the function to be differentiated. Thus, so long 
as the grid is not changed, the same differentiating macrix can be used to 
differentiate any function whose values are known at the grid points. This 
feature is shared by integrating matrices (Vakhitov, 1966; Hunter, 1970; 
Lakin, 1979) which may be used to numerically integrate functions on discrete 
grids. Differentiating matrices of the present type have previously been 
derived only for grids with equally spaced points. On these uniform grids, 
the same Newton forward difference interpolation formulas lead to both an 
integrating matrix and a differentiating matrix (Hunter & Jainchell, 1969). 
Indeed, these two types of matrices are natural complements. Integrating 
matrices have an identically zero row, and are thus intrinsically singular. 
Consequently, a differentiating matrix is necessary to undo an integration and 
acts as an approximate Inverse. 


The fact that these matrices depend only on the grid structure allows 
them to be used as approximations for differential and integral operators in 
eigenvalue problems. The result of this approximation is a matrix eigenvalue 
problem which can be solved by standard methods. 

Work in recent years has focused almost exclusively on the development 
and use of integrating matrices to solve problems associated with rotating 
beams (White & Malatino, 1975; Kvaternik, White, & Kaza, 1978; Lakin, 1982). 
Despite its links to the integrating matrix, the differentiating matrix has 
beer, almost totally ignored. Two factors contributed to this situation. 
First, for uniform grids, a differentiating matrix was found by Hunter and 
Jainchell to give less accurate results than a corresponding integrating 
matrix in a simple test problem involving the harmonic equation. Further, for 
the vibrations of a rotating cantilevered beam, the problem of principal 
interest in earlier work, it was found that the matrix eigenvalue problem 
could be formulated using only the integrating matrix (White, 1978). However, 
this is not the case for other problems of current interest. For example, in 
the rotating beam model for vibrations of spokes of an energy-storing 
flywheel, the boundary conditions are such that reformulation as a matrix 
eigenvalue problem requires the use of both differentiating and integrating 
matrices. Indeed, this is the case for many beam problems which do not 
involve cantilevered boundary conditions. 

In the present work, differentiating matrices are first derived for grids 
with arbitrarily spaced points. This is a natural generalization, as physical 
problems often have spatial variations in material properties which 
necessitate the use of unequal increments in the independent variable. 
However, in the present context, removal of the uniform grid restriction also 



serves another, more fundamental, purpose. For differentiating matrices on 
uniform grids, relatively large approximation errors are found to occur at the 
endpoints of the Interval being considered. Errors at interior points are 
considerably smaller. This situation is not surprising given the one-sided 
nature of the approximation at an endpoint and the intrinsically two-sided 
character of the derivative operation. In a formulation which is not 
restricted to uniform grids, one or more additional interior grid points can 
be added "close" to each endpoint. It will be shown that, at the cost of a 
slight increase in the size of the differentiating matrix, the addition of 
even a single "near-boundary" point dramatically reduces the errors in the 
derivative approximations. A reconsideration of the test problem of Hunter 
and Jainchell further shows that the results obtained using the 
differentiating matrix with near-boundary points as an operator in an 
eigenvalue context are now fully comparable in accuracy with the earlier 
integrating matrix results. The reliability of the generalized 
differentiating matrix is further demonstrated by considering two additional 
test problems involving both integrating and differentiating matrices. 


2. Derivation of the Differentiating Matrix 

Suppose that f(x) is a function whose values are known on a grid G 
consisting of the N+l discrete points Xq, x^,...,x N# The spacings 
xj Xj_j(j * 1,2,...,N) must be non-negative, but are otherwise arbitrary. 
Le ; {f} denote the N+l -dimensional column vector (fQ, f^,...,f^) t where 

£j * f(xj). ^ differentiating matrix on the grid G is then a square N+l- 
by-N+1 matrix [D] such that 



[D] {f} - {£'} 


( 2 . 1 ) 


where (f'J - (f', with f' - f"(Xj), is a column vector of the 

derivative values at the grid points. In particular, let [jDjJ denote the 
i-th row vector of the matrix [D], Then, by (2.1), the dot product 
|D ± J *{f) * f'. Hence, JjD^J can be obtained directly from any expression for 
f£ which can be written as a linear combination of the values fj for 
j ■ 0,...,N. In practice, the row vectors of [D] will be obtained from 
approximations to f£ rather than exact expressions. Consequently, a 
differentiating matrix on G will not be unique. 

Differentiating matrices on G may be obtained in a number of ways 
Including interpolation and least-squares polynomial approximation, and the 
form of the matrix depends on both the underlying technique and the degree of 
the approximating polynomials. In the present derivation, rows of [D] will 
be obtained using Lagrange Interpolation polynomials of degree n £ N on a 
sequence of appropriate subgrids with n+1 points. 

As a first step in obtaining a suitable approximation to fj, and hence, 
[DjJ , f(x) will be approximated by an interpolating polynomial of degree n 
on an appropriate subgrid G^ of G which contains the point x^. In 
particular, suppose that G^, consists of the n+1 grid points 
x y> x y+i » • • • ,x y+n* t * le eventua l goal is an approximation to f£, x^ 

should be centered as much as possible in G^. This is not a problem if n 
is even as x^ will be exactly centered in G if y - i - n/2. If n is 

odd, however, two choices of y are possible: y » i - (n+l)/2 or 
y ■ i - (n-l)/2. It will, of course, not be possible to center the endpoints 
xq or xjj or points x^ which either have values of y which are less than 
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zero, or values of y + n greater than N. Consequently, y must be 
restricted to the range 0 < y < N-n-1. 

If G is a uniform grid, f(x) can be approximated on by a Newton 

forward or backward difference interpolating polynomial which can then be 
rewritten as a linear combination of f^, f ^ , . . . ,f^ +n * Newton forward 

difference formulas were used by Hunter and Jainchell in their initial 
derivation of the differentiating matrix. These difference formulas may be 
generalized to arbitrary grids, if desired, using divided differences. 
However, it rapidly becomes quite difficult on non-uniform grids to perform 
the required separation of function values and grid information. Hence, the 
present derivation utilizes Lagrange interpolating polynomials which are valid 
for arbitrary grids and, as defined, separate function values and grid 
information. The required approximation to f(x) on G^ is now 

n 

fU)-£^(*> Vj (2.2) 


where the Lagrange coefficients Lj^(x) 

are given by 

| 
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Approximations to f'(x) on can be obtained directly from (2.2) by 


differentiation. Evaluation at now gives 

r i 

j-0 


(2.4) 


where 


n t } 

j i / j X , . - x , 

J * Y+j Y+m 


m=0 

m*j 


(2.5) 


I 

5 
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tyY 

In fact, the general form (2.5) of Lj 1 (x^) can be considerably simplified. 

Suppose that i * y + k with 0 _< k _< n. Then, for j not equal to k 


•i 

' l i 

f n ] 
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m-0 
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t 

U* k >j , 



(2.6) 


and for j ■ k 




(Y) 




m»0 

m#k 


x. - x . 
i Y+m 


(2.7) 


The i-th row vector |_D^J of [D] can now be obtained in a straightforward 

manner from equation (2.4). If the elements of this row vector are djj, then 

d<x - 0 for j < Y and j > N - y ~ n while d, » (x ) for 

ij i y-rm m 1 

m * 0 f • • • ) n# 


b 


A 



The introduction of near-boundary points can drastically change the form 
of the differentiating matrix on a grid. For example, if G is a uniform 
grid of five points with Xq * 0 and x^ * 4, the differentiating matrix 
based on fourth-degree polynomials is 



-25 

48 

-36 

16 

-3 




-3 

-10 

18 

-6 

1 



[D] -ij 

1 

-8 

0 

8 

-1 

• 

(2.8) 

-1 

6 

-18 

10 

3 




3 

-16 

36 

-48 

25_ 



Adding two near-boundary points 

at .01 

and 3 

.99 to 

the uniform grid 

to give a 

new ion-uniform grid of 

seven 

points. 

the 

differentiating matrix 

based on 


sixth-degree polynomials has the form 


j -102.33 102.37 

-.05 

.03 

-.02 

.25 

-.25 


J -97.69 

97.65 

.05 

.03 

.02 

-.25 

.25 


18.55 

■19.18 

-.16 J 

.12 

-.50 

6.35 

-6.18 


(D] - -8.27 

8.51 

-.89 

0 

.89 

-8.51 

8.27 

(2.9) 

6.18 

-6.35 

.50 -1 

.12 

.16 

19.18 

-18.55 


1 "* 25 

.25 

-.18 

.03 

-.05 

-97.65 

97.69 


L -25 

-.25 

.18 

.03 

.05 

-102.37 

102. 33_ 


Matrices which 

approximate higher 

derivatives 

may also 

be defined 

. For 

example, the second 

derivative matrix 

[d 2 1 

on G 

is an N+l-by-N+i 

mat rix 


such that 



[D 2 ] {f> - {f"> 


( 2 . 10 ) 







where {f"} is a column vector of values of f"(x) at the grid points. If 
n • N, then as noted by Hunter and Jalnchell, 


[D 2 ] - [D]‘ 


( 2 . 11 ) 


However, if n < N, then the second derivative matrix cannot be obtained 
simply by squaring [D] as the leading and/or trailing zeros in the rows of 
[D] lead to Inaccurate approximations. Rather, the second derivative matrix 
must be computed directly from the subgrid approximations 


f"(Xi> ■ ^ L< Y > (x^f^ 


( 2 . 12 ) 


where 


n-1 n 


-E E 


p*0 m-p+1 

P*j m*J 


2L j (Y) (x i ) 

x y+p )( Vj ■ W 


(2.13) 


As is the case with the first derivative matrix, considerable simplification 
of the coefficients in (2.13) is possible. For example, if i » y + then 


n-1 n 


’"<v -E E 


p“0 m*p+l 
p*k m#k 


Lx. - x , ;ix. 
i Y + P 1 


(2.14) 


Similar results are easily derived for third and higher derivatives 
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3. Accuracy of the Differentiating Matrix 

The accuracy of derivative approximations obtained using differentiating 
matrices based on uniform grids has been investigated by Hunter and 

Jainchell. They found that, in general, if n-th degre< polynomials are used 

to obtain approximations to the values f'(i ■ 1 N), then the errors are 

of order h n where h ■ (x^ - xq)/N is the uniform spacing. However, in 
assessing this error estimate , it must be noted that in actual applications to 
problems in rotor dynamics, h is usually 0(1) or even 0(10). In practice, 
the use of differentiating matrices on uniform grids leads to large errors in 
derivative approximations at endpoints. These errors can be clearly seen from 
the data in the first three rows of Table 1. The first row of this table 

gives the exact derivatives of f(x) ■ sin (kxx/4) and g(x) * cos(kirx/4) at 
x ■ 0 for k - 1,2,3. Rows two and three give approximations to f'(0) and 
g'(0) obtained using a differentiating matrix based on a uniform grid of five 
points with Xq ■ 0, x N ■ 4, h ■ 1, and interpolating polynomials of degrei ~ 

three and four. The errors are particularly pronounced for the two higher 

modes. 

The final eight rows of Table 1 show approximations to f'(0) and 
g'(0) obtained by using a non-uniform grid of seven points consisting of the 
previous uniform grid plus two near-boundary points at Ax and - Ax. 

Differentiating matrices based on fifth and sixth degree polynomials for Ax ■ 
.05, .01, .001, and .0001 were used. At the smallest value of Ax, 

approximations to both f'(0) and g'(0) agree with the exact values to at 
least four decimal places for all three values cl k. 

A simple example will serve to illustrate the reaaon that employing a 
non-uniform grid with near-boundary points serves to drastically reduce 



_M - l** 
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derivatlve approximation errors at the endpoints of the grid. Consider first 
the uniform grid appro 1. mat ion (with spacing h) to the derivative of an 
unspecified function f(x) at the right hand endpoint x N . Using a third 
degree Lagrange interpolating polynomial on the subgrid consisting of the 
four points Xjj _3 through x^ gives as an approximation to f' the 
expression 



1 If N " 18f N- 1 + 9f N-2 

Sh 



(3.1) 


T .f the function values fj t (k ■ N - N - 1) are now expanded in Taylor 

series relative to the endpoint x^, then 


lViJ* {£) " r( V + 0<h3) (3,2) 

which is the result found by Hunter and Jalnchell. 

Consider next the approximation to the derivative at x^ obtained using 
a fourth degree Lagrange interpolating polynomial on the non-uniform grid 

A A 

Gy consisting of G^ plus the near-boundary point x ■ - Ax. Let 

f ■ f(x). Then 


I 




By the nature of the near-bound-. ry point of the non-uniform grid, Ax is small 

A 

compared to h. Expanding the function values f and f^(lc • N-3,.,.,N-1) 
in Taylor series relative to Xjj now gives the result 


\ 

5 




- 12 - 






lVJ- (f > ■ f '<V + °(t> 7)- (3 - 4) 

Unlike error estimates for the uniform grid, the non-uniform grid error 
estimate (3*4) is essentially independent of the degree of the interpolating 
polynomials. 

An equation similar to (3.4) predicts that the errors in the derivative 

approximation at x Q will also be 0(~, A measure of the extent to 

h 

which this is valid in practice can be obtained by examining the data in Table 
1. The actual errors for the non-uniform grid approximation^ are completely 
consistent. For h * 1 and all four values of Ax, approximations to both 
f'(0) and g'(0) are valid with eirors of order Ax. 

4« The Differentiating Matrix as an Operator 

As a test of integrating and differentiating matrices, Hunter and 
Jainchell considered the harmonic equation 

y" + y * 0 0 _< x _< 4 (4.1) 

with the homogeneous boundary conditions 

y(0) * y(4) * Or (4.2) 

The exact solution of this problem is simply 

y k (x) * sin x with o) fc - . 


(4.3) 



L~. I 1 t*W 
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If equation (4.1) is written at each of the N+l points of a discrete 
grid G with Xq ■ 0 and Xj^ * 4, the differentiating matrix may be used as 
an operator to rewrite the problem in matrix notation. Let {0} be the 
N+l dimensional zero column vector. The differentiating matrix formulation 
of (4.1) is then 

U> 2 ] (y> + « 2 (y} * (0) (4.4) 

where [D 2 J is the second derivative matrix. This may be written in the 
eigenvalue form 

[A] {y} * A{y) ; [A] * [D 2 J 1 (4.5) 

where 

A * 1/d)^ (4.6) 

so that the dominant eigenvalue corresponds to the lowest frequency. Both 
boundary conditions in (4.2) have not been explicitly used in obtaining the 
eigenvalue problem (4.5). However, since y^ * = equation (4.5) reduces 

to the (N-l)-by-(N-l) dimensional problem 

(A) (7) - A{y) ; [A] « [^f 1 (4.7) 

where [D^] is the (N-l)-by-(N-l) matrix obtained from [D 2 ] by removing 
the first and last row and column and {y} is the (N-l)-dimensional column 
vector (yj, y 2 , . . . ,y N . 1 ) t . 

Hunter and Jainchell used a differentiating matrix based on fourth degree 
polynomials on a uniform five-point grid with h * 1 to solve the 3-by-3 
eigenvalue problem corresponding to (4.7). They obtained as approximations to 
the first three modes the eigenvalues and eigenvectors 
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«o - 0.7893 

l * 

»-* 

N 

O 

• 

O 

►— 

• 

o 

o 

o 

• 

oo 

S-/ 

It 

1 

| u> 2 - 1.4142 

(y 2 > - ( l.oo, o.oo, -i.oo) t . 

! w 3 - 1.7920 

| 

(y 3 ) - (-0.27, 1.00, -0.27) t . 

! 

t 

,j By contrast, exact values for 

\ 

the first three modes are 

i 

i u> 1 - 0.7854 

(yj = ( 0.71, 1.00, 0.79) t , 

‘ I 
? 

1 u> 2 * 1.5708 

{y 2 > - ( 1.00, 0.00, -1.00) t , 

5 «, > 2.3562 

.1 

{y 3 ) - (-0.71, 1.00, -0.71) t . 


The differentiating matrix results on the uniform grid thus compared poorly 
with the exact solutions* 

Hunter and Jainchell also solved this eigenvalue problem using a 
formulation baseJ only on integrating matrices on the uniform five-point 
grid* They obtained the results 

Ui l = 0.7840 {y^} * ( 0.71, 1.00, 0.71) t , 

<o 2 - 1.5492 (7 2 ) * ( 1.00, 0.00, -l.OO)*, 

co 3 * 2.8522 {y 3 > * (-0.83, 1.00, -0.83)*. 

For uniform grids, the integrating matrix formulation of the eigenvalue 
problem was thus clearly superior. 

The situation is quite different when the eigenvalue problem (4.5) is 
solved using a second derivative matrix based on a grid which includes near- 
boundary points. In particular, consider the non-uniform, seven-point grid 
obtained from the previous uniform grid by adding two near-boundary points at 


i 





I „ V 
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0.01 and 3.99. Using a differentiating matrix on this grid based on sixth 
degree interpolating polynomials, the eigenvalues and eigenvectors of the 
5-by-5 problem corresponding to (4.7) are now 


u 1 - 0.7855 
<a 2 = 1.5499 
u> 3 = 2.1724 


(y v ) * ( 0.71, 

(7 2 ) - ( i.oo, 

{y 3 } * (-0.61, 


1.00, o^dS 

0.00, -I.OO) 1 , 

1.00, -0.6 l) t . 


These results are fully comparable with the integrating matrix 
approximations. This conclusion is not changed if the integrating matrix 
method is also enhanced through the addition of near-boundary points. 

Two additional eigenvalue problems involving equation (4.1) will further 
serve to illustrate the utility of differentiating matrix operators on grids 
which include near-boundary points. Consider first the eigenvalue problem 
consisting of equation (4.1) with boundary conditions 


y'(o) = y(4) ® 0. 


(4.8) 


Exact solutions are 


/ \ A A (2k + l)ir t , 

y(x) * cos x with tu^ * g — — , (4.9) 

and thus 

Wj* 0.392/, o >2 = 1.1781, and u> 3 - 1.9635. 

Because of the form of the boundary conditions, the problem (4.1) and 
(4.8) cannot be reduced to a matrix eigenvalue problem for {y} itself using 
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only the differentiating matrix. Rather, a combination of differentiating and 
Integrating matrices can be used to obtain a matrix eigenvalue problem for the 
first derivative vector {y~}. To do this, a differentiating matrix is used 
to express the second derivative vector as {y"} “ [D] {y'}. The boundary 
condition at x^ is now used to express y(x) in terms of y'(x) as 

4 

y(x) = -/ y'(s)ds. (4.10) 

x 

If (J^] is an integrating matrix on G which approximates integrals from 
x to Xjj, (4.10) has the matrix form 

ID] (y " } - m 2 U 1 ] (y'J * {0}. (4.11) 

Equation (4.1) thus has the matrix eigenvalue formulation 

[A] {y'J - X{y'} with [A] = [D]” 1 ^] (4.12) 

A AO 

and X « 1/u . The boundary condition that y ' * 0 has not been explicity 
used in deriving (4.12). However, this condition implies solutions can be 
obtained /.am the N-by-N problem [A] (y'J ■ (y'}, where [A] is the N-by- 
N {"'•tux obtained from [A] by deleting the first row and column and 

O ) * (yj. y 2 y N > • 

To solve the matrix eigenvalue problem, G was taken to be the same 
seven-point grid as above with near boundary points at 0.01 and 3.99. A 
differentiating matrix [D] based on fourth degree polynomials and an 
integrating matrix [Jj] based on fifth degree polynomials were used. This 



even/odd degree scheme was chosen so that grid points or subintervals would be 
centered to the maximum extent possible in the approximation subgrids. 
Approximations for the first three eigenvalues were found to be 

Wj* 0.3926, o >2 * 1.1750, and <03 * 2.0883. 

Approximations to the eigenvectors {y^} for these eigenvalues, obtained by 
using the eigenvectors {y'} of (4.12) in (4.11), agreed well with exact 
results. 

In an alternate procedure, the eigenvalue problem with boundary 
conditions (4.8) can be solved for {y} directly using a combination of the 
integrating matrix [Jj] and the integrating matrix [Jq] which approximates 
integrals from xq to x. This formulation is somewhat simpler and also 
fractionally more accurate than the mixed differentiating-integrating matrix 
formulation above. However, if the boundary conditions associated with 
equation (4.1) are changed to 


y'(0) - y'(4) - 0, (4.13) 

it is no longer possible to derive a viable matrix eigenvalue formulation of 
the type tA) {y} * A{y} using integrating matrices alone. The lack of 
information on y itself at either endpoint forces use of a boundary matrix 
operator which leads to a singular matrix at the next-to-last step of the 
reformulation. This is not the case if a combination of integrating and 
differentiating matrices is used to reformulate (4.1) and (4.13) as a matrix 
eigenvalue problem for {y'}« 
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. 


-A. 


y 


If equation (4.1) Is Integrated from x to x N using the condition 
0, then 

2 4 

-y' + fa) j y(s)ds ■ 0. (4.14) 

x 


The function y(x) may also be expressed in terms of y'(x) as 


x 

y(x) = / y'(s)ds + y(0) (4.15) 

0 

and, from equation (4.1), the unknown boundary value yQ can be written as 
Yq * -Xy"(0) with X as in (4.6). Consequently, 

A s 

-y'(x) + (x - 4) y"(0) + J J / y'U)d?ds = 0. (4.16) 

x 0 

The value of y"(0) may now be obtained from y'(x) using a matrix 
formulation on a grid G involving both the differentiating matrix [D] and 
the left boundary matrix [Bq]. In particular, [Bq] is an (N+l)-by-(N+l) 
matrix with first column vector {1} and all other entries zero so that 
[Bq] {f} * {f (0) > . If [E] is a diagonal matrix with entries 
e^ * x^_j - Xjj, the matrix eigenvalue formulation of (4.1) and (4.13) is 
now [A] {y'} » X{y'} with 


[A] - {[E] [Bq] [D] + [Jq] (4.17) 


and [I] is the identity matrx. The condition that y'(0) ■ 0 implies that 
the first row and column of [A] may be removed and solutions obtained from 
the resulting N-by-N eigenvalue problem. Once and {y'} have been 
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determlned, {y^} itself may be obtained using (4.15) and (4.16) Exact 
solutions are given by (4.3) with sine replaced by cosine. 

The matrix formulation for this eigenvalue problem was solved on the 
seven-point non-uniform grid with near-boundary points G * (0., 0.001, 1., 
2., 3., 3.999, 4.). The differentiating matrix was based on sixth degree 
polynomials, and both Integrating matrices were based on fifth degree 
polynomials. Eigenvalues were found to be 

o)^ “ 0.7872, o >2 “ 1.7663, and * 3.7402. 

Results for the higher modes can be improved by increasing the number of 
points in G. 


5. Conclusions 

Generalization of the differentiating matrices to grids with non- 
uniformly spaced points removes many of the difficulties associated with 
previous versions of these matrices which were restricted to uniform grids. 
The inclusion of appropriate near-boundary points allows for accurate 
approximation of derivative values, including derivatives at endpoints. 
Greatly improved accuracy is also achieved when the present differentiating 
matrices are used as operators to express second order eigenvalue problems in 
matrix form. 

For eigenvalue problems which can be reformulated as matrix problems 
using only the integrating matrix, such as the fourth order models for 
vibrations of rotating cantilevered beams, alternate formulations involving 




the integrating and/or differentiating matrices will probably not be used in 
practice. This is because a pure integrating matrix approach, if feasible, 
allows for a slightly simpler treatment of boundary conditions. However, 
boundary conditions for many eigenvalue problems are such that a mixed 
formulation involving both types of matrices is unavoidable. Differentiating 
and integrating matrices on grids with near-boundary points should provide the 
basis for a fast and efficient numerical technique for solving these problems. 
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Table 1. Comparison of the differentiating matrices on a uniform grid of five points with h * 1 and on 
a non-uniform grid obtained by adding two near-boundary points with spacing Ax to the uniform 
grid. The test functions are f(x) * sin kxx/4 and g'x) - cos kwx/'. Approximations i~o 
f"(0) and g~(0) are obtained using interpolating polynomials of degree m as indicated. 



0.785398 0.0 1.570768 -.000012 2.356000 0.000256 



